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STABILITY OF CYLINDRICAL AND CONICAL SHELLS OF CIRCULAR 
CROSS SECTION, WITH SIMULTANEOUS ACTION OF AXIAL 
COMPRESSION AND EXTERNAL NORMAL PRESSURE* 

By Kh. M. Mushtari and A. V. Sachenkov 


We consider In this report the determination of the upper limit of 
critical loads in the case of simultaneous action of a compressive force, 
uniformly distributed over plane cross sections, and of isotropic external 
normal pressure on cylindrical or conical shells of circular cross sec- 
tion. As a starting point we use the differential equations for neutral 
equilibrium of conical shells (ref. l) which have been used for the solu- 
tion of the problem of stability of conical shells under torsion and under 
axial compression (ref. 2) j upon solution of the problem it is possible 
to satisfy all boundary conditions, in contrast to the report (ref. 3) 
where no attention is paid to the fulfillment of the boundary conditions 
and to the report (ref. 4) where only part of the boundary conditions are 
satisfied by solution of the problem according to Galerkin's method. 
Approximate formulas are used for the determination of the critical 
external normal pressure with simultaneous action of longitudinal com- 
pression. Let us note that the formulas suggested in reference 5 are 
not well founded and may lead, in a member of cases, to a substantial 
mistake in the magnitude of the critical load. 


1 . SYMBOLS 


We shall use the following symbols: 

2y angle of taper 

r distance, along the generatrix, between the vertex of the 

cone and a point on the median surface 

Tq distance to the smaller one of the bases 


* 
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L 

2h 

9 

n 


w 


Po 


length of the shell along the generatrix 

thickness of the shell 

angle "between the axial plane and the plane under 
consideration - 

number of waves forming along the periphery when the shell 
"buckles 

normal displacement of the point on the median surface 
during buckling 

external norma], pressure, acting both on the lateral surface 
and on the bases of the shell section considered 

additional compressive force applied to the smaller of the 
end cross sections 


m m 
10 * 20 


T r T 2 , S 

e 2 

E 

cr 


linear membrane forces up to buckling, determined according 
to the momentless theory 

additional membrane forces appearing during buckling 
additional annular expansion 
modulus of elasticity 
lateral-expansion coefficient 


K = 


2Eh 
1 - cr 2 


D = 


2Eh- 


3(l - a 2 ) 


(rigidities of the shell) 


l (1.1) 


r = 1 + t = In t m-i = (m- integer) 

r 0 z 


Rq = r 0 tan 7 


V - (1 - °> 

2 


H-L = 1 - 2v 2 

„ - n 
1 sin 7 


2 4 

= v - 3v 


9^ = cp sin 7 


( 1 . 2 ) 
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€ = 


Vh 2 (l - t 2v - 2 ) 


3 r Q 2 cot 2 7(1 - a 2 )^ - v) (t 2v - 1 ) 


D' = 


(l + 2v)Eh(T 2v - l) 


A = 


jvr Q ta y( T ( 1+2v ^ - l)J 

(1.3) 

T n (l + 2v) (l - t^ 2v_1 ^) 0.5 + A ^ 2 + 0.25(1 + 2v) 2 ] 

X 1 = 


PqRqCI - 2 v)(t( 1+2v ) - 1 ) 

e 2 (m 1 2 + v 2 ):^ 






Tl 2 = 


m-L 2 + (1 - v) 2 J + H}^ 2 + 


K* = 


D' 

n ^ 2 + 0.25(1 + 

2 v) 2 ] 

( mi 4 

+ Him! 2 +• ) 



(m 2. 2 + v 2 ) 

] 




(1.4) 


2. FULFILLMENT OF THE BOUNDARY CONDITIONS AND INTEGRATION 
OF THE EQUATION OF COMBINED DEFORMATION 


tt-ip differential equations after introduction of the stress function 
will he (ref. 1 ) 


Dr£Aw - cot 7 ^ - r 


dff 

dr 2 


m d^W , m /l Sw + _L 
%0 + T 20(r 3r r2 i(f 2l 


= 0 


( 2 . 1 ) 


,2 

r *AAf + 2Eh COt 7 = 0 

Sr 2 
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# 


where 


A( ) . afLL + i MJL + i 

Sr 2 r 3r r 2 ^2 


T _ = iM+I 5 2 f 

r Sr r 2 Scp-^ 2 


To = 


Sr 2 


S = 


. dfj 

3r\r 3cf^j 


( 2 . 2 ) 


_ p 0 r tan 7 Vo T _ „ r tan r 

T 10 5 1 20 P0 r tan y 


(2.3) 


If the ends of the shell are simply supported, the conditions 


v =* 0 


. a(bw , 1 i „ q 

3r 2 r V* r 3cp 2 


(2A) 


= 0 - 0 


must be satisfied on the edges r = rg and r = Tq + L. 

We Introduce the following substitutions: 

z = In ~ f = F cos n^cp^ v = eVZw- 3 _ cos n-^cp^ t = In 

(2.5) ^ 

Then, taking equations (2.3) into consideration, we bring equations (2.1) 
into the form 



2Ehr 0 cot ye 


V. 
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5 (v-4)zl 


d\r 


4z^ 


^ - 4(l - v) 


dw 


^ 3 . + K] + Kp^ ^ + K,_v, - 

dz3 ^ dz 2 ^ ^ 5n 1 


r 0 COt l a (l-v) Z , 


/ d 2 F dr\ 

Vdz 2 ‘ **) 


PO tajl 7 , 5z 

+ — 5 V e ^ 


,2 

i iii 

2 dz 2 


( v+ 5 


\dw 

i/dz 


W eZ 

D 

where 


+ ( 2 „ - i) + (v2 - v) J 


^ = 0 


K la = 6v 2 - 12v + 4 - 2n 1 2 

K 2 n = 4v 3 - 12v 2 + 8v + 4n 1 2 (l - v) 


K 


3n 


= v^ 1 " - 4v ^ + 4v 2 - n^ 2 ^ - 4v + 2v 2 ) + n-j^" 


~ _ 2 l 1 2 

-n i " 2 2 V 


K 4n w l I 


( 2 . 6 ) 


(2.7) 


The "boundary conditions (2.4) assume the form 

d 2 ^ 


w-j_ = 0 

- n T 2 ! 1 = 0 

dz - 1 - 


= 0 

dz* 

d 2 F _ dF = 0 
dz 2 


For z = 0 and z = t 


( 2 . 8 ) 


(2.9) 


We shall seek the solution of the boundary- value problem, assuming the 
shape of the wave formation 


Wi = Aq sin m-^z 


( 2 . 10 ) 


Thereby, the boundary conditions (2.8) will be satisfied 
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From the first equation of the system (2.6) we find 


F = A^e 11 -^ 2 + Age n ^ Z + B^e 


(ni+2)z + B 2 e (2 - n l )z + 


2A^EhTQ cot /e ^ +V ) Z Am 0 1 B m l 2 + Xnm C0B ra l z ) (2.11) 


where A l> 


B^, Bg are arbitrary constants 




ran 


\in 


%tn 

^mn 


( m l 2 + v - v2 )%m “ 4 < v ’ 2 ^)% 2 *ran] : (%in 2 + l 6 m ±^m?) 

*1 j(l - 2v)% n + kvfef + v - v 2 )*,^ : (q^ 2 + l6m 1 2 v 2 ^ mn 2 j 

(m x 2 + ni 2 ) 2 -^2(1 + v 2 )^ 2 + 2(1 - Jv^m-L 2 + (v 2 - l ) 2 

m^ 2 + nj 2 + 1 — v 2 

( 2 . 12 ) 


Satisfying the conditions (2.9) , we arrive at the system of equa- 
tions for determining Aj_, Ag, B^, and Bg. 


As has been shown in reference 1, the initial equations assumed were 
used only in the case of a thin shell of small length, the buckling of 
which occurs with the formation of a considerable number of waves along 
the circumference. 


For the approximate determination of the hyperbolic terms in the 
expression (2.11) we can take 

sinh n-[_t cosh n^t 


especially since, as will be shown later, those terms affect the magnitude 
of the critical load only insignificantly. In this manner, we find 


*6 + n lS 


A 1 « . _ u e ( L+v )‘ fc 


2 ^nj_ 2 - n^) sinh n^t 




B l « IS e ( v - 1 ) t 

2(n± + l) sinh n^t 


B2 “ 


% - n i $ 5 


n i + n i 


2 

L 

*5 


> (2.13) 


n ! - 1 
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where 


*5 - 0 . 2 5 (^ - nfxj) 

$lj. =* (1 + v) ®g + m-|_$-j_ 

*1 “ C 1 •+ v ) V. “ m l X mn 


*6 - °-5( 2 *2 * "A* * \) 

= (l + v)$q_ - 

% = C 1 + + m lVi 


(2.14) 


3 . THE CHARACTERISTIC EQUATION 


Substituting equations (2.10) and (2.11) into the second equation of 
system ( 2 . 6 ), we integrate that equation according to the method of 

Bubnow-OalerkLn, multiplying it by e( 2+v ) z sin m-jh Az- 


in this manner, we arrive at the characteristic equation 


P-Q,-R+T = 0 


(5.1) 


where 


P = 


Q - 


R = 


- m l g 2 [ m l*~ + ^( v ~~ 1 ) 2 ]a l 2 ~ m l^ln + ( 1 ~ v ) K 2 n + K 3 n] 

m-^ 2 + (l - v 2 ) 

- 1. 5v(l + v) - 0.25 + n^ 2 ^l” 1 ! 2 + 3v(l ~ v ) ~ O. 5 J 


Ptf 11 ! 


4v 


(•t 2v - l) 


m -, 2 + 0 . 25(1 + 2 v) 
2 vt 


m x 2 + 0 . 25(1 - 2 v ) 2 


(*6 * n l°5) e V ( 4 6 ~ n i $ 5) 

m l 2 + ( n l + v - 1 f m i 2 + (n x + v - l) 


»5 e2Vt (*l + 2 ) T *5^1 - 2 ) 


+ (n^ + 1 + v ) 2 m-j 2 + - 1 - vy 


T = m iK - %) + v K • $ 2) 

m -, 2 + v 2 


(5-2) 



8 


NACA m 1433 


This equation may "be greatly simplified for the wide range of thin 
shells which satisfy the condition 


t - 



L < 1.72r 0 


Then we have according to equations (1.2), m-i 2 ^ n 2 . In addition, 

O " L 

for thin shells of small length n^ is much greater than 1. Thereby 

we have, retaining only the main terms in the expressions of the quanti- 
ties (2.12), (2.14), and (3.2) 


kv 

on^ - Zon-jp - m-j_2 

n^ 2 + 2(l - a)m^ 


2 

m l 

m-|_ 2 + (n-^ + v + l) 2 J 

m-j 2 +“(112 + v - l) 2 


Calculations according to this formula, which take" into account the 
expression P and the solution obtained further on, show that the quan- 
tity R in equation (3»l) may be neglected if we permit an error in the 
magnitude of the critical load of 2 to 3 percent in the„ direction of 
increased load. In the remaining terms of the equation we can also 
neglect the quantities of the order one in comparison with n-j_2. This 
gives an additional error of 1 to 2 percent. 


Let us note that the maximum error is admitted on the boundary of 
the region where the shell becomes rather long and, for stability, 
requires transverse reinforcements. With the reduction of t 1116 per- 
missible error rapidly drops to 1 to 2 percent. Thus, we shall determine 
the critical load from the approximate equation 

P q ( ftp 2 + 0.5mi 2 fo 1 ! 2 \ _ 

D ' \ m l 2 + 0.25(1 + 2v) 2 m x 2 + 0.25(1 - 2v ) 2 j 

e 2 (m^ 2 + nq 2 ) 2 mp^ + Pimp 2 + P2 
m-L 2 + (l - v) 2 (m-j 2 + n 1 2 ) 2 ( m 1 2 + v 2 ) 

In the special case of longitudinal compression 
than 1 and taking the designations (1.3) and (1*4) for Pq = 0 into 
consideration, we find 


T 0 (l + 2v)< 

(l - t^-D) 

_ € 2 I 

( m l 2 + n l 2 ] 

I 2 

_ 1 _ . 

mi 2 

D'(l - 2 v ) 1 

( T (^v) . 

m x 2 ! 

[-1 2 + H! 2 ) 2 


(3-3) 

m^ 2 is much greater 
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Consequently, the critical compressive force is equal to 

T 0 - (1 - g -> ^ -/y_ ( 5 .4) 

Jv(l - v) Rq(i - T^ v "^) 

The deviation of Tq given by the above formula from the corre- 
sponding exact formula obtained by Shtaerman does not exceed 4 to 5 per- 
cent for t = 1. 


4. DETERMINATION OF THE CRITICAL LOAD 
Using the notation (1.4) and, in addition 


6 = 


2 2 
m-1 + nq* 

mi 


we bring equation (3*5) into the form 

From the condition dpg^db = 0 we find 


2 ft 4 

ti 5 = 


, + 2 - jH 

• 

1 + 2 ( A - l) — ^ 

3 + 5 

• 

1 + ^Aq lj g 


(4.1) 


(4.2) 


(4.3) 


By means of simple, but rather tedious calculations we can show 
that dpg/dmq is greater than 0. 

Such a monotonic increase of the quantity Pq with increasing m^_ 
can be explained by the fact that the quantities K', Aq, and are 

changed very little by the increase of m^, owing to the fact that, as 
can be seen from equation (4.2), for every fixed & the minimum Pq is 
reached at the largest of the values permitted by the boundary 

conditions . Therefore 


(4.4) 


m = 1 


m-L = n/t 
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From = 1 in equation (4.3) we find 

d = Qj_ = \B : fi\ 


(4.5) 


o 

If » 1, this solution holds in the case where the prehuckling 

stress is equal to the ring stress. 

In the general case 


5 « 


Si 

1 + P 


(4.6) 


where 3 satisfies the equation 

2(1 - _ [l - (1 + P)^| 

8 l (1 + 3)[l - i(l + 3) 4 " 


Hence it follows that 


-0.2 ^ 3 £ 0.2, if 0.86 > 2^1 - > -2.93 


(4.7) 


(4.8) 


Thus we may assume for 3 the smaller (with regard to the modulus) 
root of the equation 

e a [s> - M 1 - >i)§] + n[« - 2 4 - + 2 (! - - 0 0 *- 9 ) 

and the approximate value of the critical pressure is equal to 

P 0fk = 1.31K'Tl 3 / 2 m 1 (l.33 + 23 2 ) : 1 + ^(\ - l)(l + P)j (4.10) 


In the neighborhood of its minimum, the value Pq changes slowly; 

therefore, the critical value of the pressure determined for 3 (which 
was found from (4.9)) differs even on the boundaries of the region (4.8) 
from its value for 3 which satisfies the condition of minimizing (4.7) 
by less than 0.6 percent, although the error in the magnitude of 3 
attains 13 percent. “ 
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For shells which satisfy the condition 


0.49 £ 2(1 - X x )|i £ -0.83 


(4.11) 


that is, for -0.1 ^ (3 ^0.1, admitting a conservative error of about 
1 percent, we may put (3 = 0 in equation (4.10) . In this manner, we 
arrive at the simple formula 


P0,k 


1.74K'm 1 T) 5 / 2 

1 + (mi/6i)(\i - 1) 


(4.12) 


Here and in equation (4.10), the quantities K', m-|_, 7}, &n<i 

are determined according to formulas (4.4), (4. 5 ) , and (l.l) to (1.4). 

For X = 0, X x = 0.5 we obtain therefore the following formula 
for the case of isotropic pressure 


®6,k" 


1.74K'm 1 T] 5 / 2 
1 - 0.5m-j_/6^ 


(4.13) 


In the special case of a cylindrical shell with the radius R we 

have 


7 = 0 t « k , sin 7 _ _ 

R 1 


nR 

L sin 7 


K' = 


2Eh 


R sin 2 7 


X, - 0.5 + X X = — 


_ h Bin^y 


P 0 R R\p(l - a 2 ) 


According to formula (4.12) we find the magnitude of the critical 
isotropic pressure in the case of simultaneous action of the axial com- 
pressive force Tq: 


4.8gEh / h' 
L(l - o 2 ) 5 /4\R 


3/2 


1 + 1.82 /hi . / T 0 

L(1 - a 2 ) ^ V P 0 R 



PO = 


(4.14) 
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whereby the condition 


O .265 > (1 - 2 A)ip S’ -0.45 

Lt 


(4.15) 


must be fulfilled. 

Let us note that equation (4.14) is equivalent to the equation 


where 

m _ 2.68Eh 2 

°“ "e/oT^ 

This last equation is approximately satisfied only under condi- 
tion (4.15). Moreover, it is necessary to take into consideration the 
case where the shell may lose its stability "with a bang" due to axial 
compression. Therefore, our formulas ought to be used only when the 

total axial force Tq + ^PqR is smaller than the lower limit of the 
critical axial force T -| H which, as is known (ref. 6), is determined 
according to the formula 

TlH ~ ( 4 . 18 ) 

Let, furthermore , Tq + ^PqR be less than or equal to 0. 78 Eh 2 /]R 
Then, according to equation (4.l6) 

PO = POmf 1 " 


!-&*2- 

P 0 R 


0.78Eh 2 ' 


HP, 


1 - 


Om 


p 0m R \ 
^Gm ) 


(4.19) 


1.56Eh g (l - Pp^/gp 
P0nI R 2 (l - 0.78Eh 2 /m 0m ) 


T 0,m p 0,m 


(4.16) 


_ _ 4.85Eh (h 

' L(1 - ^5/W 


5/2 


1 _ 0 . 91 /hR 

L(l - a 2 ) 1 ^ 


(4.17) 


Consequently, 
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or, taking equations (4.17) into consideration, assuming a = 0.3 and 
neglecting in the second-order terms h R/L^ in comparison with unity, 
ve have 

(1 - g? — - 0.42 > -0.45 

Ii Xj 


that is, condition (4.15) is fulfilled. 


Formal a (4.12) which determines the critical normal pressure in 
the case of simultaneous action of axial compression, can "be greatly 
simplified under the condition ^ n. For 0.25 ^ a ^ 0.53 (as is 

usually the case with metals), taking into account new notation, 
formula (4.12) after easy calculations and neglect of second-order 
terms, can he put in the form 

= , g (2 - <r) - l) 1 ^! - T -11 *' 1 ) ^ E tan 7(h/RpA 

' /IT7(l - c^KrC 2 - a > - l)ln t P + »(* - 0-5)] 

(4.20) 


where 


T _ 1 + L_ 

^ (2 - ff)T 0 (l - r" a ) 

r 0 

PqV 1 

(r(2-a) . 1) 

8 - ^ . 

'[l - " 

lA 

(~h tan 7 

+ O 

- l] | 

^Rq Ua T 


(4.21) 


Let us note that in the change from formula (4.12) to the simplified 
formula (4.20) we reduced the critical pressure at most hy 5 percent 
(for t = l); hut formula (4.12), in turn, was derived from the character- 
istic equation (3*l), hy way of simplifications, which increased the 
pressure we had heen seeking hy 2 to 5 percent. In this manner, formula 
(4.20) ultimately gives the upper limit of the critical pressure, reduced 
hy less than 3 percent. 


Moreover, we must not forget that formula (4.12) and notahly, also, 
formula (4.20) are derived for the range of variable shell parameters 
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determined, by the inequalities (4.11) which, after simplification, 
assume the form: 


0.49 ^ (1 - 2A)e £ -0.83 


(4.22) 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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